On the Casimir entropy for a ball in front of a plane 



M. Borda^l 

Leipzig University, Institute for Theoretical Physics, 04109 Leipzig, Germany 

I.G. PirozhenkcQ 

Bogoliubov Laboratory of Theoretical Physics, 
Joint Institute for Nuclear Research and Dubna International University 
Dubna 141 980, Russia 
(Dated: October 2, 2010) 

The violation of the third law of thermodynamics for metals described by the Drude model and 
for dielectrics with finite dc conductivity is one of the most interesting problems in the field of the 
Casimir effect. It manifests itself as a non-vanishing of the entropy for vanishing temperature. We 
review the relevant calculations for plane surfaces and calculate the corresponding contributions for 
a ball in front of a plane. In this geometry, these appear in much the same way as for parallel planes. 
We conclude that the violation of the 3rd law is not related to the infinite size of the planes. 

PACS numbers: 03.70.+k Theory of quantized fields 
ll.10.Wx Finite-temperature field theory 
11. 80. La Multiple scattering 
12.20.Ds Specific calculations 



a bordag@itp.uni-leipzig.dc 
b pirozhen@theor.jinr.ru 



2 



I. INTRODUCTION 

The violation of the third law of thermodynamics by the Casimir effect for certain properties of the interacting 
bodies is still one of the most interesting problems in the field. It manifests itself as a non vanishing limit of the 
entropy S, defined as minus the derivative of the free energy T with respect to the temperature T, 

for T — >• 0. It was first observed in [l|, 0] for metals described by the Drude model and in Q for dielectrics with dc 
conductivity. Inserting the corresponding permittivities into the Lifshitz formula for the free energy J 7 , at vanishing 
temperature a term linear in T remains which by means of ([1]) results in a non vanishing contribution to S at T = 0. 
Obviously, the use of these permittivities, which otherwise work fine, when inserted into the Lifshitz formula, results in 
a behavior which is not only non acceptable from a principle point of view but which is also in quite clear disagreement 
with experimental observations. Details on this topic can be found in the recent review Q (and also in the book [f|). 
The Drude model is characterized by a permittivity 

^) = i + ^4tT' (2) 

where uj = i£ is the imaginary frequency, uj p is the plasma frequency and 7 is the relaxation parameter. For 7 = 0, 
the permittivity e D (i£) turns into that of the plasma model which does not cause problems with thermodynamics. A 
dielectric with dc conductivity is characterized by a permittivity 

e dc (iO = co(»0 + , (3) 

where cr is the dc conductivity and £o(*£) is the permittivity of a dielectric without dc conductivity of which we only 
need to know that it takes a finite limit eo = eo(0) for zero frequency. 

In both cases the violation occurs if the parameters cto or 7 are non zero, depend on the temperature and decrease 
for T — ¥ 0. This happens for some reasonable idealizations of real materials where <Jo decreases exponentially fast or 
7 as a power of T (for metals with perfect crystal lattice). 

In the present paper we reconsider the derivation of the violation terms and establish that an arbitrary slow decrease 
already results in a violation of the third law. For this purpose, we employ a representation different from that used 
in H which looks technically more direct and which allows for a deeper insight into the corresponding structures. We 
would like to stress that we make no claim about the physical reality of slowly decreasing parameters 7 and u. 

During the past decade there was quite a number of attempts to avoid a violation of the third law. Most of them 
point to a modification by including additional physical effects. An example is the addition of impurities to a perfect 
crystal lattice Q ■ Other consist in using impedance boundary conditions in place of the Drude model in the Lifshitz 
formula It must be admitted that no satisfactory understanding was reached so far. 

In the present paper we answer the question whether a finite size of one of the interacting bodies is able to prevent 
the violation. For this we consider a sphere in front of a plane at low temperature. This is an extension of our 
previous paper [Uj . We consider a sphere with the permittivities fl2J and ([3|) in front of a conducting plane. Special 
attention is paid to the case of small separation. We mention that the configuration of a ball in front of a plane at 
finite temperature is under active discussion, see for example [12[ and [13[. 

In the next section we reconsider the case of parallel planes and in the third section we treat the sphere-plane 
configuration. Conclusions are drawn in the last section. 
Throughout the paper we use units with h = c = 1. 

II. THE FREE ENERGY FOR PARALLEL PLANES 

The free energy for two parallel plane bodies at separation a is given by the Lifshitz formula 

^ E i E Ml-r^e—) (4) 



2 ^ J (27r) 

i=-oo v ' i=TE,TM 



with 



yJ# + &, (5) 
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the Matsubara frequencies 

6 = 2tt/T, (6) 
and k is the momentum parallel to the planes. The reflection coefficients are 

for the TE mode and 



r eg- y (£ -l)ff + g 2 
eg + v/( e - 1)£ 2 + g 2 

for the TM mode, where for e one must insert one of the two, ([2]) or ([3]), according to the model considered. 
Using the Abel-Plana formula, representation (j4]) can be rewritten as a sum, 

T = E + A T J", (9) 

of the vacuum energy, 

J ° i=TE.TM 



depending on T only through 7 or a, and the explicitely temperature dependent part of the free energy. 

1 

where 



A T J"=-^/ dxn T (x) $(a:), (11) 



is the Boltzmann factor. In (|lll) the notation 

$(x) = i((p(ix) — (p(—ix)) (13) 

is used which contains the two contributions appearing from turning the integration path £ — ► zbix in the Abel-Plana 
formula and 



p(0 = / dkk J2 k(l-r?e- 2o ») 1 (14) 

•'° i=TE,TM 

is the function which must be analytically continued. This representation of the free energy is well known and it is 
especially useful at low temperature where the Matsubara sum in (0J) converges slowly. 

It must be mentioned that the division in @ is done according to the contributions of the electromagnetic exci- 
tations. The vacuum energy Eq does not contain contributions from the thermal photons whereas AtF is just their 
contribution. However, in the case of temperature dependent permittivity, the thermal excitation of the degrees of 
freedom of the interacting bodies like electronic or phonon excitations enter the vacuum energy through the cor- 
responding parameters "f{T) or cr(T). In this way their vacuum energy contributes also to the entropy which thus 
consists of two parts, 

S = S + Si (15) 

with 

_ dE dE dp 
S °-—9T -—gfaf' (16) 



where p is one of the two, 7 or a, and 



9 - dATT (U) 
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The first part, Sq depends on the thermal excitations of the interacting bodies only, whereas the second part, Si, 
depends on the thermal excitations of the photons too. 

The more conventional approach to derive the violation terms, used in the literature (see the book [5| for a de- 
tailed representation), rests on the observation that the interesting terms result from the (I = 0)-contribution to the 
Matsubara sum in ((4]). As it turned out, for metallic bodies described by the Drude model the linear term is 



T 

j--— = -j-j~>- + ... =-- 

with 



-7-Drudc -j-plasma.TE . rB f-t o\ 

T = ~^=° + - (18) 



/£, = (2a) 2 ^ dqq In (l - (r™ ma ) 2 e -*«) . (19) 

Here rj^, m is the reflection coefficient of the TE mode with the permittivity of the plasma model, i.e., e D (i£), ©, 
with 7 = 0. For large uj p the linear in T contribution takes the limiting value 



^ Dmd V p _ = T^C(3) + .... (20) 
lD7rcr 



For a dielectric body with dc conductivity, the contribution linear in T can be written as a difference between the 
(I = 0)-contributions of the TM mode with dc conductivity and without (<tq — in Eq.Q), 

x-dc _ x-dc.TM t-iio dc,TM . _ T Dc 

* -• / "/=0 -^1=0 +■■■ ~ ~ 167] . a2 /lin +■■■ ( ZS ~) 

with 

/,? n C =C(3)-Li 3 (r 2 ) (22) 

and the notation 

-trl < 23 » 

for the TM reflection coefficient for static permittivity. Li3 is a polylogarithm function. 

It is interesting to remember that the (I = 0)-contribution to the Matsubara sum gives just the leading order 
contribution for T — > oo. In this way, the violating terms are closely related to the high temperature limit, i.e., to the 
classical limit. It is, however, not clear whether this has any deeper meaning. 

In the following subsections we use representation © with (|10p and (fTTj) and re-calculate the low temperature 
behavior in both models, Drude and dc conductivity. 

A. The vacuum energy in the Drude model 

In this subsection we calculate the contribution of the vacuum energy, E (-y), to the entropy at small temperature 
which comes in through the temperature dependence of the relaxation parameter 7(T). According to (|16p we have 

s °- — d^Tdf- (24) 

We assume a decrease of 7(T) for T — > according to 

7 (T) = 7l T« + . . . (25) 

with a > 0. Therefore we need the expansion of the vacuum energy for small 7. As it follows from the calculations 
below, this expansion contains a logarithmic contribution, 

£0(7) = £ (0) + 7 (- ln(2a 7 ) Ex+E^+.... (26) 

Here, E o (0) is the vacuum energy of the plasma model and it does not contribute to S^. The logarithmic term results 
from the TE mode. With ([TO]!. © and © we note 



£o TE (7) = i r f # K (?) , (27) 



71TE 

1-./,, -"Jo '"' Vi" 
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where we defined 



h q (z) = In 



1 - 



-2aq 



(28) 



As compared with (|10[) we changed the integration over k for q = t; 2 + k 2 and interchanged the orders of the 
integrations. 

In (|27p . a formal expansion for small 7 goes in powers of 7/^ and already in the first order the ^-integration becomes 
singular. For this reason we change the variable £ for z = 7/f, 



-1 />oo />oo 



and integrate by parts two times 
^o TE (7) 



1 

4^2 



dqq < qh q 
V9 



-*".iMi 



7 / dz \nzh q (z) 



(29) 



(30) 



Here the primes denote derivatives of the function h q (z) with respect to z and we used the property of this function 
and of its derivatives to decrease sufficiently fast for z — > 00. In representation (l30l) . it is possible to expand up to the 
first order in 7 and we get for the expansion parameters in (|26[) 



E 1 = ^L 1 ,/,,,/,;„„ 



and 



E r = dqq 



(l+\n(2aq))h'J0) 



dz lnz h„ (z) 



(31) 



(32) 



In both expressions, the integrations are well convergent and deliver certain smooth functions of the plasma frequency 
uip. We show Ei as a function of ui p in Fig. Q] (left panel). For large u> p it decreases, 



£1 



C(3) 



O 



1 



(33) 



It should be mentioned that this logarithmic contribution was already found in eq.(17). 

The contribution from the TM mode is easier since in that case it is possible to expand the integrand directly into 
powers of 7 not producing singularities in the integrations (at least in first order in 7). As a consequence, there is no 
logarithmic contribution and Ej M is a smooth function of uj p like Ej E and it can be calculated easily numerically. In 
this way, in fllJJ} we have to insert E 1, (EU, and E 1 = Ej E + E™. 
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FIG. 1. The coefficient E\ in (|26|l for the Drude model, (|3ip . as a function of lo v for a separation a — 1 (left panel) and for a 
dielectric with dc conductivity, (|45|l . as a function of eo (right panel). 
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B. The vacuum energy for a dielectric with dc conductivity 

The calculation of the vacuum energy for a dielectric with dc conductivity goes in principle similar to that in the 
preceding subsection. In place of (|24[) we have now 



■Si 



dE (a) da 
da~df' 



(34) 



which depends on the temperature via the conductivity cr(T) in the permittivity ([3]) (we dropped the factor 47T and 
defined a = 47rc7o) and we assume 



a(T) = aiT a 



(35) 



with a > 0. Again we will observe a logarithmic contribution for T — > 0. This time it comes from the TM mode and 
we note a formula in parallel to ([26)1 . 



E (a) = Eb(Q) + a (- ln(2acr) E x + E x 



(36) 

for small a. Starting from here, for technical reasons, we proceed in a slightly different way. We start from the TM 
contribution and consider its derivative, 



d 



l 



with the notation 



da 



^™ = 4^ / dqq / =±g[?,t,q) , 



i - 



(eo +z)q- V(eo - 1 + z)£ 2 + <? 2 



-2aq 



In (|37|) we integrate by parts in the variable £, 



_5 
d 



) -E™(a) = 1 ^ J™ dqq |ln(2o«) 3 (j, 9, <7 



,f>9 



(37) 



(38) 



(39) 



where we used that g(a/l;,£,q) vanishes for £ — > 0. 

The integral in the figure brackets in this equation requires a special trea tment. We divide it into two parts, A\ 
and A2, according to the division of the integration region into £ G [0, \faja\ and £ € \y/o~fa, q\. This is possible for 
any fixed c; since we are interested in the limit a — > 0. So we have 



A x = 



I '' a «M2«o|*(f,e,*) 



and A2 accordingly. In this integral we change the integration variable for £ = £/<r, 



Ax 



d( ln(2acrC) ^sIt,^,? 



Now it is possible to tend cr — > and we get 



At 



d( ln(2ac<) 



f 

d( ° vc 



:,0,g + 



(40) 



(41) 



(42) 



where the dots denote higher powers in a. In the contribution proportional to the logarithm of cr, the integration can 
be carried out using the derivative, 



A, =]n(2aa)g(0,0,q) + J d(ln( JLg(l,0,q 



(43) 



with 



7(0,0,?) 



(44) 
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and ro is the reflection coefficient (|23|) . Since this is the only place where a logarithm of a appears we can now write 
down Ex in (ffi)]) . With (1591 and the above two formulas the integration over q can be done and we obtain 

- _ 1 U 2 (r 2 ) 

El -^{l-el)a?> (45) 

where Li 2 is the poly logarithm. Next we need to calculate A 2 . This is quite easy since we can put a = directly in 
the integrand, 

A 2 = [\hx(2a0^z g(P,t,q) • (46) 
Jo d£ 

Collecting all contributions we get from ([39]) and (|37|) the contribution from the TM mode to E\ , 

-i poo f />oo a / -i \ 

E™ = —J o dqq hn(2aq)g(0 >q ,q)- J d( ln( — g I -, 0, q) 

- jf° d£ ln(2a# ^ <?(0, £,<?)} . (47) 

The integrals in this expression are convergent and deliver a smooth function of eo which can be evaluated numerically. 
The contribution from the TE mode does not have a logarithmic term and it can be calculated by simply putting 
a = in a formula in parallel to (|57|) . The emerging integrals are finite too. We restrict ourselves here with a 
representation of E\ as a function of eo for a = 1 in Fig. [T] (right panel). 

In this way, we observe in both models a logarithmic contribution to the vacuum energy for small parameter, 7 or 
a. If we insert these, together with (|25p resp. ([55")) , into (IM1) rcsp. we obtain, say for the Drude model, 

Su = -a ll T a ~ 1 ({a\nT + l)E 1 +E^j + ... , (48) 

and a similar formula for dc conductivity. We observe not only a non zero, but even a diverging contribution to the 
entropy at vanishing temperature for < a < 1. 



C. The temperature dependent part of the free energy in the Drude model 

We start the calculation of the temperature dependent part of the free energy from eqs. (|TT|) and (fl3|) with the 
function <p(£), eq. ([13]). with £ = ix inserted. So we have to consider 

/>oo 

tp(ix)= dkk Yj hi (1 - r? e- 2o «) , (49) 

~' i=TE,TM 



with g = vP-?- We divide the integration into a first region, fc g [0, x], and a second region, k £ [x, 00). We need 
the function tp(ix) for small x ~ T. Therefore, in the first part of the integration region we have k < T and a factor 
~ T 2 from gJ/c/c. As a consequence, the contribution from this region is by two additional powers of T suppressed as 
compared with the second region where the integration region is infinite. Hence a contribution to the linear in T term 
can come from the second region only. In that region we change the variable of integration from k for q = \Jk 2 — x 1 
(which is real) and arrive at a representation of the temperature dependent part of the free energy, up to higher orders 
in T, given by eqs. (TrTj) . ([121) and (TT51) with a function 

/>oo 

<p(ix)= dqq ln(l-rf e- 20 *). (50) 

"'O i=TE,TM 

The reflection coefficients are still given by eqs. ([7J and © but with the substitution —> ix and with the permit- 
tivities ([!]) and with £ — > ix. 

Now we consider the Drude model. For small T, the leading contribution results from the TE mode. With ([2]) and 
10) we get 

iT ^, TE = _l /; _^ /; dqqi (,„ $ _ ,„ (3.)) + . . . (51) 
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with the function h q (y) defined in (|28|) . Now we consider a > 1 in (|25|) . i.e., 7 decreases not slower than the first 
power of the temperature. We make the substitution x = 7C 



A jrDrude.TE 



1 

4^2 



7 



e ii/T _ x 



dg g z ( h. 



-'X 



and note the ratio 7/T in the Boltzmann factor. Now we can tend T 
we get 



in the integrand and note (p25)l . For a = f 



^ jrDrude.TE 



167ra* 



■/ D (7i) 



with 



/ D ( 7 i) = (2a) 5 



d(- 



7i 



' eTi< - 1 

This is a smooth function of 71 and of 2au> p . For a > 1, because of 



'/'/'/' j /',, ( — ) - /'., ( 



(52) 



(53) 



(54) 



for T — > 0, we get from 
Its explicit expression reads 



the linear in T term (|5^|) with the same function / D , however with zero arg ument, / D (0) 



n»> = ff 



dqq 



XI 



In 1 - 



c.c. 



(55) 



where we inserted the explicit expression ([28|i and substituted g = g/2a. Here, and below, we use the notation 'c.c' 
indicating that the complex conjugate inside the square bracket must be subtracted. We would like to comment on 
the know property of the linear in T term not to depend on the relaxation parameter. This property holds for a > 1 
and it results from the ratio in Eq. (|54[) having a smooth limit for T — ¥ which does not depend on 71. 

In the other case of a decreasing relaxation coefficient where < a < 1 holds, we start from ea. (l5"lT) and make the 
substitutions x 



= TC and q = q^/T/j 
T 2 



A T J" 



Drudc.TE 



47T 2 7 



d(- 



7 



7 



7 



(56) 



In this representation it is possible to tend T — > in the integrand. The emerging integrals are finite. However, the 
factor in front, T 2 /-f ~ T 2 ~ a goes to zero faster than a first power of T. Hence it does not contribute to the violation 
of the 3rd law. 

It is easy to see that / D (0), Eq. ([55"T) . is a smooth function of 2auj p . It is also possible to show that this is just the 
same function as /° n , Eq. (fH))) in ([15]): / D (0) = /j? . We have plotted / D (0) as a function of cjp (with a = 1) in Fig. 
H (left panel). 

In this way, for a decrease of the relaxation parameter j(T) linear in T or faster we have a linear contribution to 
the free energy. If the decrease of 7(T) is slower this linear term disappears. Finally we remind that the TM mode 
does not contribute to a linear term. 

For a better understanding of the structures involved and as reference for the case of a sphere in front of the plane 
we remind here shortly the case of a fixed relaxation parameter 7. In that case there is no linear term and T 2 is 
the leading order which receives contributions from both, TE and TM modes. The starting point is again eq. (|51j) 
with the function (z) = h q (z), (|28[). for the TE mode. The contribution from the TM mode is given by the same 
formula, ([5"Tj). with 
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instead. For the TE mode one needs to make the substitution x = TC, and q = yTq. After that one can put T = 
directly in the integrands and one comes to the expression 



A jrDrude,TE 



T2_ 

4tt 2 



sC-1 



dqq i < In 



1 - 



+ <7 2 



(58) 



These integrals can be computed easily delivering 



A x"Drudo,TE = 2 In 2 - 1 cj 2 T 2 
T 48 7 



(59) 



We note that this expression is just the same as that which follows from (|53|) for large 71 with the formal substitution 
71 -> 7/T. 

For the TM mode one starts from the same formula, (|510 , now with (|57p inserted. Here one needs only to substitute 
x = TC, and one can expand the integrand for small T. This expansion starts from a first order term which delivers 



A jrDrude,TM 



T 



4tt 2 J eC - 1 J 



dqq 



-8-fxT 
(e 2a <? - 1) uj 2 



Carrying out the integration one comes to 



A j7Drudc,TM _ 



7T 2 7 T 2 

" 18 (2a) 2 w 2 



(60) 



(61) 



which must be added to (159 



D. The temperature dependent part of the free energy for a dielectric with dc conductivity 



For a dielectric with dc conductivity we apply the same scheme of calculations as in the preceding subsection. In 
this case the interesting contribution comes from the TM mode and in place of (|5l"j) . using ([8|) and <j3j> , we have now 

At-F™ = ± jf" [dq ql [Xn(l- r 2 TM ^) - ex.] (62) 

with 



(eo(iar) + 5) 9 ~ Vpeo(^) + 1) x 2 + ixcr + q 2 
(e {ix) + §) g + ^(-eo(ia;) + l).T 2 +ia;(7 + Q 2 

We make the substitution x = <rC an d consider cr(T) decreasing according to (pJ5j) with a > 1. After this substitution 
we can tend T — > in the integrand. For a = 1 we come to 

Ar7 -dc,TM = _ T fd 

lb7ra z 
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with 



/ d >i) = - 



(2a) 2 

4tt 2 



0"! 



:<7iC - 1 



dqqi 



In 1 




— 2aq 



This is a smooth function of o\ taking a finite limiting value for u\ = and it decreases for a\ 
come to the same expression (|64p . but with the function f dc with zero argument, 



(65) 

oo. For a > 1 we 



/ dc (0) = - 



1 

4^2 



dqqi 



In 1 - 



1 + zC(eo - 1) 
1 + i({e + 1) 



In 1 



1 - iCjeo - 1) 
1 - zC(e + 1) 



(66) 



So we see in this case a picture similar to the Drude model. If the conductivity a(T) decreases at least linear with 
T, we have a linear in T contribution to the temperature dependent part of the free energy. If it decreases slower 
(that is, a < 1) this linear term can be shown to be absent. 

In (|66|) . in the reflection coefficient, the g-dependence dropped out and we substituted q = q/2a. This expression 
can be simplified by an expansion of the logarithm into a sum after which the q-intcgration can be done immediately, 



r c (o) 



_L_V— f°° dx . 



fc — 1 



1 + ix(e — 1) 
1 + ix(e + 1) 



2k 



1 - ix(e - 1) 
1 — ix(e + 1) 



2// 



Now the x-integration can be done too and we arrive at 

/ dc (0)=Li 3 (r 2 )-C(3). 



(67) 



(68) 



This is just the same function j22J as in eq. (pT|) . i.e., / DC (0) = holds. We show it on Fig. ([2} (right panel). 
Since in this case the TE mode does not contribute we reproduced the linear term (|T5)) . 

We conclude this subsection by comparing with the case of a fixed DC conductivity a. The temperature dependent 
part of the free energy is given by eq. (|5ip with the functions 



h™(a,x) =ln 



eo-l + SU 2 



(69) 



and 



h™(a,x) = In 



1 - 



'(^o + fh-^o-i + D^ + g 2 



-laq 



X £ o + f)9+V(^o-l + f 



r 



(70) 



for the two modes. In both cases one needs to make the substitution x = TQ. In the TE case one needs to integrate 
by parts for q several times to see that the TE mode results in a order T contribution and can be dropped. In the 
TM contribution one can expand 



h™ 



8CT 



cr(l - e 2a i] 

and carry out the remaining integrations. The result is simply 



+ 



(71) 



TM 



T 2 



18 (2a) 2 a 



(72) 



which is the counterpart of (fBT]) . Moreover, ([72^1 and (fBTj) are related by the substitution a — > j/uj 2 which just relates 
e D , $Z§, with e dc , ([3]), in leading order for £ —> (note our notation a = 47ro- )- 

It should be mentioned that the low temperature behavior for fixed parameters in known. It was calculated earlier 
using other methods. For example, Eq. (|5"9"]) can be found in [l4| and Eq. (1721 in [l5j |. 



11 



III. THE FREE ENERGY FOR A SPHERE IN FRONT OF A PLANE 

In this section we extend the results of the preceding section to the geometry of a sphere, metallic or dielectric, in 
front of a conducting plane. The radius of the sphere is R and and the distance from its center to the plane is L. 
Technically, this is an extension of our previous papers, [ll[ and [l6j]. So for the temperature dependent part of the 
free energy we use the same basic formula, 

1 f°° 

A T T = — dxn T (x) iTr[ln(l-M(ix)) - In (1 - M(-ix))] , (73) 



2^ Jo 

with the Boltzmann factor given by Eq. (fT2)) . The matrix M(£) has the entries 

with the coefficients H l w and A\ v given by eqs.(7) and (10) in [ll[ and it is a matrix with respect to the polarizations. 
For the TE mode the function dj B (£) is given by 

where n = y/ejl is the refraction index. The function dJ M (£) for the TM mode follows by interchanging e and fi. 

These functions are expressed in terms of the know modified spherical Bessel functions, si(z) — y/ 7rz/2/; +1 / 2 (.z) 
and ei(z) = y/2z/irKi +1 / 2 {z). For the following it is useful to separate the powers of the argument in front of the 
ascending scries. We represent 



si(z) = V^r(§) +1 *z(z), ei(z) = Vtt (f)' +1 fc/(z), 



(76) 



with 

»,(*) = (f)"^ 1/2 h{z) = |(|) i+1/2 ^ +1/2 (z), 

il(z) = + l + »,(*), ki{z) = (-l + z&)ki(z). 
All these functions have power series expansions, e.g., 

ii(z) = H (0)+i i l 1) z 2 + ... , ~k(z) = i l (0) + i< i 1) z 2 + ... , (78) 



(77) 



which we will need below 



These definitions allow to rewrite the functions dj (£) in the form 



2 / < \ 2l + 1 



with 



and 



= jgi«(Oi«(ygo-i«(Qi«(^) 

where we restricted ourselves the a pure dielectric ball putting ji = 1. Now the entries of the matrix M take the form 

^ (eR/2) w +1 / A r A \ / t TB (efl) o \ , . 
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We redistributed a factor (£R/2y which is admissible under the trace in ((73)) . 

For the low temperature expansion this expression can be simplified. Below we need the lowest orders for £ — > 0, 



W + ... • (83) 



From the powers of £ in ({82} it follows that only I" = I + I 1 contributes to (|82|) and that Mo is diagonal in the 



polarizations. The latter follows form the additional factor of £ in Aj For this reason, and because of the trace in 



(|73p. the temperature dependent part of the free energy becomes a sum of the two polarizations, 

AJ - = AJ 7 ™ + AJ 7 ™ + . . . , (84) 

which holds in the orders of T we are interested in. We continue with a separate consideration of the Drude model 
and the dc conductivity in the following two subsections. 



A. Ball described by the Drude model 

In this subsection we consider a ball described by the Drude model in front of a conducting plane. It turns out that 
we can act in quite close analogy to the planar case in subsection 2.3. We assume a relaxation parameter decreasing 
with temperature according to (|25|) . We make the substitutions x — 7C in (|73|) and with ([2]) we note 



1 and 



e D (7C) 



\A D (7C) 7C = ^y^+--. (85) 



for 7 — » 0. After that we expand the matrix M for small 7. First we expand i TE (£). With (|85l) we get from (|80|) for 
the TE mode 



i TE (70=*o E (C,^) + ... (86) 



with 



*0 (C^pj- —/ r— s 1 7 r— \ 1 I 8 ') 

whereas for the TM mode 

i™( 7 C) = 0( 7 2 ) (88) 

holds. So, only the TE mode contributes to the linear in T term. In the remaining factors in M we can put also 7 = 
and come to 

M^ udc (C)^M M ,( 7 0| T .o, 

= X (2 ) 2 ' +1 ^(O^A^^WpiD , (89) 

where we defined 

P=j (90) 

for the ratio of radius of the ball to the distance between the plane and the center of the ball. It must be mentioned 
that M°™ do (C) still depends on ( through tJ E because of (g5J). 
Now we insert M ; D ™ dc (C), ©, into (J73]), 



d( —^jl iTr [(1 - M Drudc (iC)) - (1 - M Dmdc (-<))] . (91) 



This formula is in parallel to (|52j) in the planar case. It has a non vanishing limit for T — ¥ in case of a 7 decreasing 
not slower than linear, i.e., for a > 1 in (|2"5)l . Now we restrict ourselves to a > 1 and come to 



At J" = f - f&afaupR) + 0(T 2 ) (92) 
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FIG. 3. The function /bailCPi w p)> (|93|) . as a function of p for several values of uj p (left panel) and the function /baii(Pi e o), (|99 
as a function of p for several values of eo (right panel). 



with 



/ b D all (p, Up R) = J- [°° ^ £Tr [(1 - M D ™ dc «)) - (1 - 



trDrude 



HO)] 



(93) 



This function describes the contribution linear in T in the Drude model. In Fig|3] we have plotted f^aiiiP'^p) as 
function of p for several values of u p . For the calculation of the trace one needs to make a truncation of the orbital 
momenta, I < l m . In this case it turned out that for all values of the parameters p and io p a few lowest I < 4 are 
sufficient. This includes, for instance, the case p = 1, i.e., zero separation. 



B. Dielectric ball with DC conductivity 

In this subsection we derive the linear contribution for a dielectric ball with DC conductivity. We use the permittivity 
(J3J) and assume a decrease of the conductivity according to (1331) . We substitute £ = cr£ and note 



1 



c 



, r/r 



K) + C 



and Je dc ((T() a^R = <7^/((( + e )R + .... 



(94) 



We have to insert these expansions into tf E (£), (l80|) . and £™(£)> <[5T|>. For the TE case we observe that i™(crC) is of 
higher order in er and, for this reason, it does not contribute to the linear in T term. In opposite, for the TM mode 
we have 



*o m (^)^ m W) M 



^^(0)1,(0) -*,(0)i,(0) 



(95) 



which is the analog of eq. (|87p for the Drude model. The remaining calculations go in parallel to the preceding 
subsection. In the remaing parts of the matrix M we put a — and come to 



MtAO = ^ (f) 2i+1 fc ;+ ,(0)^'A|+^™( C ,a) 
Now we insert this Mf j, (£) into (J75J) and get 

1 



A T J" = 



- i-: 



(7c 



HO)] 



(96) 



(97) 



In this way, also for dc conductivity we have a linear contribution for T — > in case a vanishes linear in T. If it 
vanishes faster, just like before, the linear term becomes independent form a and it is 



A T F=^-f^(p,e ) + O(T 2 ) 



(98) 
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with 

ftUP, 60) = J™ J *Tr [(1 - M rfc «)) - (1 - M dc HC))] . (99) 

This function can be calculated in the same way as f^ a n(p, uj p ) in the preceding subsection. We have plotted f}Zn(p, eo) 
in Fig. [3] as function of p for several values of e . We used a truncation I < l m of the orbital momenta and again some 
low I were sufficient. 

C. Both models with fixed parameters 

In this subsection we consider both models, Drude model and DC conductivity, with fixed parameters. The starting 
point is again Eg. (1751) for the temperature dependent part of the free energy. For the limit T — > 0, we need the matrix 
M(£) for small £. First we consider the permittivities for the Drude model. With @ and we note 

-4 = ^ + ..., V^R=^v / ?+---, (ioo) 

for £ — > 0. For the DC conductivity using ^ we have 

-L = -£ + •••> V^SR = y/ZRy/i+..., (101) 

We see that in this approximation both models are related by the substitution a — > Wp/7. Therefor we can restrict 
ourselves to the Drude model. 

We need to inserted ()100|) into (|80|) and (|81|) . In lowest order in £ we get 

™=*r^+... (102) 



with 



7 



.TE wr-wi 



for the TE mode, and 



with 



tr = " ' * y * , (io3) 



to 1M (0=*o 1M +*r^ + --- (104) 



TM *i(0) / inc \ 



*z(0) ' 



.tm _ *i(0) f 1 , *i(0)*i(0) 



t M = ^ -1 + :' v ; 1V ; , (106) 

ki(o) V fci(o)*i(o), 



and 



for the TM mode. We insert these expressions into the matrix M, (J82J), separately for both modes. For TE mode 
there is no zeroth order and up to the first order we find 

uj 2 R 2 

M TE = M™ — i + ... (107) 

7 

with 



M 1 T V = ^(f) W ' +1 ^K0)^'A|/*F. (l«).s) 
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For the TM mode we have 



,TM 



1™ 



,TM 7 



(109) 



with 



M™ 

JW 1,1' 



2 V2/ '-•-r» \~'~ 1,1' 1,1' 

t (f) ,+, ' +1 ^(o)«/r< t ™. 

We have to insert these expansions into the logarithm in (|T3[) . Expanding the logarithm we get 



u> 2 R 2 

Trln (1 - M TE (0) = -TrM^ E £ - 

Finally we insert this expression into (j73[) and carry out the ^-integration, 



(110) 



(111) 



uj£R 2 T 2 



(112) 



with 



(> 



TrM™, 



(113) 



which is the leading order in the low temperature expansion of the TE mode contribution to the free energy in the 
Drude model with fixed parameter 7. The function g TE (p) can be calculated numerically. For that one needs to 
truncate the orbital momenta I < l m . The emerging expression turns out to be converging for l m — > 00 for all 
p € [0, 1]. The function g TE (p) is shown in Fig. @] (left panel). In fact, Eq. (|113p gives a power series expansion of this 
function. This can be seen from (|108[) and the absence of M in this case. 
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-1.0' 




FIG. 4. The function <? TE (p) for several values of the truncation parameter, l m — 1, ... ,11 (left panel) and the function g™(p) 
for several values of the truncation parameter, l m = 1, . . . , 8 (right panel). The dashed lines corresponds to l m = 1. The limit 
of small separation corresponds to p = 1. 

The corresponding expression for the TM mode has to account for a nonzero M™ given by Eq. (|110[) and the 
expansion for £ — > reads 



Trln (1 - M™(£)) = Trln (l - M^ M ) - Tr (l - 



ttTM\ 



(iTM 



Being inserted into (|T3")) only the odd term survives and after the ^-integration we get 



A T T™=g^(p) 



TM/ 



IT 2 



(114) 



(115) 



with 



TM 



(P) 



•Tr (1 



I™) 



irTM 



(116) 
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This function can also be calculated numerically making a truncation as before. However, in this case the convergence 
is weaker. More exactly, for any fixed p < 1 there is convergence for l m — » oo, for small p > it is even very fast 
converging. For p = 1 there is no convergence. The function g TE (l) grows not slower than linear with l m . We 
interpret this as a non-commutativity of the limits T — > and p — > 1. As a consequence we expect for p = 1, i.e., 
for contact, a slower decrease with T as in (|115[) . However, the case of contact is unphysical since in that case the 
vacuum energy is infinite. We have plotted the function g™(p) in Fig. |4] (right panel). 



IV. CONCLUSIONS 



In the preceding section we have shown that for a ball in front of a plane there is a violation of Nernst's theorem 
(3rd law of thermodynamics) in much the same manner as for parallel planes. For the Drude model, with a relaxation 
parameter 7 decreasing faster than linear with the temperature T, the entropy at T = is given by 

s u = ~ fZi(P^ P R) (H7) 

Z7T 

with the function f^ ail (p, ui p R), Eq. (|93|) . shown in Fig. [3] (left panel). For a dielectric ball with dc conductivity a also 
decreasing not slower than linear with T, the residual entropy is 

S dc = -^/baii(P^o) (118) 

Z7T 

with the function f^\\(p, £o), Eq. (|99|) . shown in Fig. [3] (right panel). These residual entropies are in complete analogy 
to the planar case for which the corresponding functions are shown in Fig. [2] From here our conclusion is that the 
violation of the 3rd law is not related to the infinite extend of the parallel planes. It must be mentioned that this 
conclusion, strictly speaking, does not apply to the Drude model since in a finite size body the relaxation parameter, 
which is inverse proportional to the electronic mean free path, does not vanish. 

In order to derive (|117|) and (|118|) we used the representation ([73| for the temperature dependent part of the free 
energy as it appears from applying the Abel-Plana formula to the corresponding Matsubara sum. On that way we 
re-derived the violating terms for the planar case reconfirming and simplifying the original derivations, Q and 0. 
At once this allowed to extend the parameter range for which the violation occurs. We remind the two sources of 
entropy, Sq and <Si, as defined in Eq. (|15[) . Sq appears from the dependence of the vacuum energy on the temperature 
through the relaxation parameter j(T) or the conductivity a(T) and Si comes from A^T which is the temperature 
dependent part of the free energy for temperature independent 7 and a. We found that Si has a violating term if 7 
or a decrease for T — > not slower than linear, the cases 7 ~ T and a ~ T included. The contribution to the entropy 
Sq is present if the decrease is not faster than linear, i.e., 7 ~ T a or a ~ T a with < a < 1. In fact, So even diverges 
for T — > 0. We do not discuss the question whether this has any physical relevance. Our point is simply to show what 
happens if inserting such parameters into the Lifshitz formula or its generalization to more complicated geometry. 
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TABLE I. The sign of the entropy and the contributing modes for all cases considered in this paper. 
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We remind that the behavior of the entropy is completely different if 7 or a have a finite limit for T — ¥ 0. In that 
case the temperature dependent part of the free energy is proportional to T 2 and the entropy vanishes in the limit as 
it should. This is well known for parallel planes and we showed here that it holds also for a ball in front of a plane. 

As observed already several times, the Casimir entropy may take negative values. In the Table 1 we collect the cases 
considered in this paper and show the sign of the entropy for T close to zero. In addition we show the responsible 
mode(s) and the relevant formulas and figures. In all cases considered, the TE mode gives a negative contribution to 
the energy and the TM mode gives a positive one. In most cases only one mode contributes (the other goes with a 
higher power of T), in some cases both modes contribute to the leading behavior for T — > 0. In that cases both signs 
are possible independence on the parameters involved. 
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